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ABSTRACT 

A  dynamic  initialization  scheme  for  tropical  cyclone  structure  and  intensity  in  numerical  prediction  systems 
is  described  and  tested.  The  procedure  involves  the  removal  of  the  analyzed  vortex  and,  then,  insertion  of 
a  new  vortex  that  is  dynamically  initialized  to  the  observed  surface  pressure  into  the  numerical  model  initial 
conditions.  This  new  vortex  has  the  potential  to  be  more  balanced,  and  to  have  a  more  realistic  boundary  layer 
structure  than  by  adding  synthetic  data  in  the  data  assimilation  procedure  to  initialize  the  tropical  cyclone  in 
a  model.  The  dynamic  initialization  scheme  was  tested  on  multiple  tropical  cyclones  during  2008  and  2009  in 
the  North  Atlantic  and  western  North  Pacific  Ocean  basins  using  the  Naval  Research  Laboratory’s  tropical 
cyclone  version  of  the  Coupled  Ocean-Atmosphere  Mesoscale  Prediction  System  (COAMPS-TC).  The  use 
of  this  initialization  procedure  yielded  significant  improvements  in  intensity  forecasts,  with  no  degradation  in 
track  performance.  Mean  absolute  errors  in  the  maximum  sustained  surface  wind  were  reduced  by  approx¬ 
imately  5  kt  for  all  lead  times  up  to  72  h. 


1.  Introduction 

In  operational  atmospheric  numerical  prediction  sys¬ 
tems,  there  are  many  factors  that  contribute  to  forecast 
error — or  the  mean  absolute  value  of  the  difference 
between  the  model’s  prediction  of  some  quantity  (e.g., 
winds  or  pressure)  at  some  lead  time  and  the  observed 
value  of  that  quantity  at  that  future  time.  The  causes  of 
errors  may  come  from  four  main  areas  and  combinations 
of  them:  (i)  imperfect  initial  conditions,  (ii)  insufficient 
model  resolution,  (iii)  limits  of  the  representation  of 
physical  processes,  and  (iv)  limits  of  predictability.  Since 
short-term  numerical  weather  prediction  is  an  initial 
value  problem,  accurate  specification  of  the  initial  con¬ 
ditions  is  essential  to  reducing  forecast  errors.  However, 
there  will  always  be  errors  in  the  initial  conditions  be¬ 
cause  the  atmosphere  is  never  perfectly  observed,  es¬ 
pecially  over  open  oceans. 
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At  present,  there  has  been  a  significant  effort  made 
to  improve  tropical  cyclone  (hereafter  TC)  intensity 
forecasts  in  numerical  prediction  systems.  Intensity,  as 
defined  here  is  the  National  Hurricane  Center  (NHC) 
criterion  of  the  maximum  1-min  sustained  wind  at  10  m 
above  the  surface  in  the  hurricane,  or  the  minimum  cen¬ 
tral  surface  pressure.  Hurricane  intensity  forecast  skill  has 
shown  little  improvement  over  the  past  20  years.  For  ex¬ 
ample,  the  official  NHC  forecast  error  at  a  48-h  lead  time 
was  approximately  15  kt  (or  7.7  m  s-1,  where  1  kt  = 
0.514  ms-1)  for  2009  in  the  North  Atlantic  basin,  only 
a  modest  improvement  from  the  1990  value  of  16  kt  (in¬ 
formation  online  at  http://www.nhc.noaa.gov/verification/ 
verify5.shtml;  see  the  least  squares  trend  line).  During  this 
same  time  period,  TC  track  errors  have  been  reduced 
steadily.  In  2009  the  official  track  error  at  a  48-h  lead  time 
was  approximately  75  nautical  miles  (n  mi;  or  138.9  km, 
where  1  n  mi  =  1852  m)  in  the  North  Atlantic,  which  is 
a  significant  improvement  from  the  1990  value  of  ap¬ 
proximately  200  n  mi.  One  reason  for  the  lack  of  intensity 
forecast  skill  improvement  is  that  TC  intensity  is  governed 
both  by  its  interaction  with  the  large-scale  environment  as 
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well  as  internal  dynamical  processes  occurring  inside  the 
storm.  Recent  research  has  shown  that  very  high  hori¬ 
zontal  grid  resolution  (order  of  200  m-2  km)  is  needed  to 
adequately  model  some  of  these  internal  dynamical  pro¬ 
cesses  to  obtain  better  intensity  forecast  skill  (Davis  et  al. 
2008;  Gentry  and  Lackmann  2010).  With  the  limitation  of 
current  computer  resources,  it  is  generally  not  feasible  to 
run  operational  TC  models  at  such  high  horizontal  reso¬ 
lution,  and  therefore  it  cannot  be  expected  that  these 
models  will  be  able  to  fully  capture  internal  processes 
partly  responsible  for  the  intensity  change.  Nonetheless, 
there  is  room  for  improvement  in  numerical  prediction 
systems  at  these  coarser  resolutions,  and  one  area  that  we 
target  here  is  obtaining  a  better  set  of  initial  conditions. 

Due  to  the  paucity  of  observations  in  the  inner  cores 
of  most  TCs,  it  is  generally  accepted  that  vortex  bogus- 
sing  is  needed  to  improve  the  representation  of  the  TC 
in  numerical  models  (Leslie  and  Holland  1995).  There 
are  a  number  of  different  methods  that  have  been  pro¬ 
posed  for  initializing  TCs  in  numerical  prediction  systems. 
One  method  is  the  use  of  variational  data  assimilation 
with  synthetic  observations  of  a  TC  vortex  that  closely 
matches  the  observed  TC  intensity  and  structure  (Serrano 
and  Unden  1994;  Zou  and  Xiao  2000;  Xia  Pu  and  Braun 
2001;  Xiao  et  al.  2006;  Wu  et  al.  2006;  Liou  and  Sashegyi 
2011).  However,  at  present  the  balance  constraint  in  three- 
dimensional  variational  data  assimilation  (3D VAR)  sys¬ 
tems  is  often  mostly  geostrophic,  and  thus  they  are  generally 
not  able  to  produce  a  gradient-balanced  vortex.  This  can 
manifest  itself  in  the  production  a  TC  with  too  weak  a 
central  pressure  for  the  winds,  or  vice  versa,  causing  an 
initial  balance  adjustment  in  the  early  stage  of  the  nu¬ 
merical  integration.  Another  method  entails  inserting 
an  idealized  balanced  vortex  into  the  model  with  a  wind 
field  that  matches  the  observations  (Mathur  1991;  Leslie 
and  Holland  1995).  A  major  limitation  in  this  approach 
is  that  the  boundary  and  outflow  layers  in  TCs  are  not 
balanced,  so  there  will  be  an  initial  adjustment  after  the 
integration  starts.  The  method  we  advocate  here  is  a  dy¬ 
namic  initialization  method,  where  Newtonian  relaxation 
terms  are  added  to  the  right-hand  side  of  a  given  prog¬ 
nostic  variable  (e.g.,  pressure  or  velocity)  equation,  and 
a  full-physics  model  dynamically  initialized  to  the  desired 
final  intensity  in  idealized  conditions.  This  procedure  has 
the  potential  to  produce  a  more  dynamically  and  ther¬ 
modynamically  consistent  vortex,  with  a  more  realistic 
boundary  layer  and  outflow  layer,  and  the  model  is  less 
likely  to  go  through  a  period  of  shock  when  such  a  vortex 
is  initialized  in  it.  Early  studies  (Hoke  and  Anthes  1976; 
Hoke  and  Anthes  1977)  demonstrated  the  potential 
utility  of  using  a  dynamic  initialization  technique  in 
numerical  prediction  systems.  Kurihara  et  al.  (1993)  and 
Bender  et  al.  (1993)  described  a  dynamic  initialization 


scheme  using  an  axisymmetric  model  and  demonstrated 
its  improved  performance  in  the  Geophysical  Fluid 
Dynamics  Laboratory  (GFDL)  hurricane  model.  Peng 
et  al.  (1993)  demonstrated  the  utility  of  a  dynamic  ini¬ 
tialization  scheme  in  the  limited  area  primitive  equation 
model  used  by  Taiwan’s  Central  Weather  Bureau. 

The  purpose  of  the  present  work  is  to  evaluate  a  dy¬ 
namic  initialization  scheme  used  in  conjunction  with  the 
3D  VAR  system  in  the  Naval  Research  Laboratory’s  me- 
soscale  TC  prediction  model.  A  companion  study  (Zhang 
et  al.  2011,  manuscript  submitted  to  Wea.  Forecasting , 
hereafter  ZLGPP)  examines  the  dynamic  initialization 
scheme  with  the  Weather  Research  and  Forecasting 
Model  (WRF).  The  outline  of  this  paper  is  as  follows.  In 
section  2,  the  mesoscale  TC  prediction  model  is  described, 
and  both  the  control  and  dynamic  initialization  procedures 
are  discussed.  In  section  3,  structure  and  intensity  compar¬ 
isons  between  the  control  and  dynamical  initialization 
methods  are  first  presented  for  a  few  case  studies  and  then 
intensity  and  track  errors  are  provided  for  a  large  sample 
of  cases.  A  summary  of  the  study  is  presented  in  section  4. 

2.  Forecast  model  and  initialization  methods 

a.  Mesoscale  TC  prediction  model 

The  mesoscale  model  used  here  is  the  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System.  A  description 
of  the  original  CO  AMPS  model1  is  provided  by  Hodur 
(1997)  and  more  details  can  also  be  found  in  Chen  et  al. 
(2003).  The  model  uses  a  terrain-following  sigma-height 
coordinate  and  the  nonhydrostatic  compressible  equations 
of  motion  (Klemp  and  Wilhelmson  1978).  The  micro¬ 
physics  scheme  is  based  on  Rutledge  and  Hobbs  (1983), 
with  prognostic  equations  for  mixing  ratios  of  cloud 
droplets,  ice  particles,  rain,  snow,  graupel,  and  drizzle. 
The  model  also  includes  a  shortwave  and  longwave  ra¬ 
diation  scheme  (Harshvardhan  et  al.  1987),  and  a  plan¬ 
etary  boundary  layer  scheme  with  a  1.5 -order  turbulence 
closure  (Mellor  and  Yamada  1982). 

The  tropical  cyclone  prediction  version  COAMPS-TC 
includes  the  following  enhancements:  (i)  synthetic  wind 
and  mass  observations  of  the  TC  based  on  the  opera¬ 
tional  warning  message,  (ii)  relocation  of  the  first-guess 
field  to  the  observed  TC  position,  (iii)  TC-following  nested 
inner  grids  using  an  automatic  TC  tracker,  (iv)  dissipative 
heating  (Jin  et  al.  2007),  and  (v)  a  surface  drag  coefficient 
that  approaches  2.5  X  10-3  for  wind  speeds  exceeding 
35  m  s_1  (Donelan  et  al.  2004).  The  forecast  system  also 
has  a  capability  for  ocean  coupling;  however,  for  this 


1  CO  AMPS  is  a  registered  trademark  of  the  Naval  Research 
Laboratory. 


652 


WEATHER  AND  FORECASTING 


Volume  26 


36  °N 


24  °N 


12  °N 


Fig.  1.  The  COAMPS-TC  grid  setup  in  the  western  North  Pacific  basin.  The  outermost 
domain  has  a  horizontal  resolution  of  45  km  and  is  held  fixed,  and  the  two  inner  domains  (with 
horizontal  resolutions  of  15  and  5  km,  respectively)  move  with  the  TC. 


study  the  model  was  run  in  stand-alone  atmosphere 
mode. 

For  the  experiments  conducted  here,  three  nested  grids 
were  used,  with  horizontal  resolutions  of  45, 15,  and  5  km, 
respectively.  The  15-  and  5-km  grids  automatically  move 
following  the  TC  circulation  with  two-way  interactive 
nesting,  while  the  45-km  grid  was  held  fixed.  The  setup  of 
all  three  domains  for  COAMPS-TC  in  the  western  North 
Pacific  basin  is  shown  in  Fig.  1.  The  model  was  run  with  40 
sigma  levels  in  the  vertical  with  a  top  at  31  km,  and  an 
update  cycle  of  12  h  was  used.  While  the  microphysics 
scheme  was  used  for  all  three  nests,  the  Kain-Fritsch 
cumulus  scheme  was  activated  in  the  45-  and  15-km  do¬ 
mains  to  help  resolve  subgrid-scale  convection. 

For  the  first  forecast  of  a  TC,  the  global  analysis  from 
the  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS;  Hogan  and  Rosmond  1991)  is  used 
as  the  first  guess,  defined  here  as  a  cold  start.  Subsequent 
runs  of  the  same  storm  then  use  the  previous  COAMPS- 
TC  forecast  as  the  first  guess,  identified  here  as  warm 
starts.  In  either  case,  a  3DVAR  system  (Daley  and  Barker 
2001)  is  used  to  optimally  blend  the  observations  with 
the  first  guess,  creating  the  analysis.  For  TC  prediction, 
the  first  forecast  is  generally  a  cold  start,  while  all  sub¬ 
sequent  forecasts  of  it  will  be  warm  starts. 

b.  Control  initialization  procedure 

In  the  control  initialization  method,  the  TC  structure 
and  intensity  are  included  in  the  analysis  using  synthetic 


observations  based  on  the  TC  warning  message.  An 
axisymmetric  vortex  is  constructed  with  a  modified 
Rankine  wind  profile  that  fits  the  observed  TC  intensity 
and  size  parameters  of  the  maximum  wind,  and  the  radii 
of  34-  and  50-kt  winds.  From  these  parameters,  the  ra¬ 
dius  of  maximum  wind  and  the  modified  Rankine  vor¬ 
tex  decay  parameter  are  uniquely  determined  (Liou  and 
Sashegyi  2011).  Next,  the  geopotential  and  temperature 
fields  are  determined  by  enforcing  gradient  and  hydro¬ 
static  balance  constraints  on  the  winds.  Zonal  and  me¬ 
ridional  velocities,  geopotential  height,  and  temperature 
synthetic  observations  are  created  at  eight  azimuthal 
segments,  and  at  radii  of  0.5°,  1°,  2°,  4°,  and  6°.2  The 
synthetic  observations  are  included  at  the  vertical  levels 
of  1000, 965, 925, 850, 775, 700, 600, 500,  and  400  hPa,  with 
prescribed  vertical  decay  factors  of  1.000,  0.996,  0.992, 
0.983,  0.970,  0.950,  0.920,  0.870,  and  0.720,  respectively. 
Finally,  the  3DVAR  system  is  used  to  optimally  blend 
the  TC  synthetic  observations  and  all  the  other  obser¬ 
vations  with  the  first  guess  to  obtain  the  initial  fields  of 
the  model.  In  the  3VDAR  system,  the  geostrophic  bal¬ 
ance  constraint  is  relaxed  in  the  tropics  to  help  produce 
a  more  gradient-balanced  vortex.  More  details  on  the 
control  initialization  procedure  can  be  found  in  Liou 
and  Sashegyi  (2011). 


2  For  TCs  with  a  maximum  sustained  wind  of  less  than  23  m  s  1, 
the  6°  synthetic  observations  are  not  included. 
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c.  Dynamic  initialization  procedure 

The  tropical  cyclone  dynamic  initialization  procedure 
(hereafter  TCDI)  involves  the  removal  of  the  analyzed 
vortex  from  the  initial  field  of  the  numerical  model,  and 
then  the  insertion  of  a  new  vortex  that  is  dynamically 
initialized  to  the  observed  surface  pressure.  A  flowchart 
describing  how  the  TCDI  system  is  applied  to  such  a  nu¬ 
merical  prediction  system  is  shown  in  Fig.  2.  Note  that 
here  the  dynamic  initialization  scheme  is  invoked  after 
3DVAR;  however,  this  scheme  may  also  be  used  before 
3D  VAR  to  improve  the  first  guess  prior  to  the  analysis 
(ZLGPP).  The  removal  of  the  analyzed  vortex  is  done  by 
using  a  Tukey  window  spatial  low-pass  filter.  The  filter 
cutoff  wavelength  was  set  to  5  X  105 * * * * 10  m  for  zonal  and 
meridional  velocity  fields,  and  3  X  105  m  for  the  tem¬ 
perature,  geopotential  height,  and  moisture  fields. 

For  the  specific  experiments  here,  we  use  the  Tropical 
Cyclone  Model  (TCM3;  Wang  2001)  to  spin  up  the  vor¬ 
tices  in  idealized  environmental  conditions.  This  model 
has  been  previously  shown  to  produce  realistic  TC  vor¬ 
tices,  capturing  many  aspects  of  the  observed  structure 
(Wang  2002a, b).  The  model  was  run  in  hydrostatic  mode 
and  uses  a  terrain-following  sigma-pressure  coordinate.  An 
initially  weak  vortex  is  first  specified  in  the  model,  and  then 
the  vortex  is  nudged  during  integration  to  the  observed 
surface  pressure  at  the  center  of  the  model  domain  as 

5  =  ~r(X  “  ^Obs)’  (!) 

where  y  is  the  relaxation  coefficient,  ps  is  the  prognostic 

model  surface  pressure,  and pobs  is  the  estimated  surface 
pressure  from  the  best-track  data.  It  can  be  readily  seen 

that  in  the  absence  of  other  forcing  to  the  surface  pressure, 
ps  will  exponentially  decay  to  pobs  with  a  He  damping  time 

of  1/y.  For  the  experiments  here,  y  was  set  to  4.4  X 

10-5  s"1,  corresponding  to  a  1/e  damping  time  of  6.25  h. 
This  rather  small  coefficient  is  used  in  order  to  gently  relax 
the  vortex  to  the  prescribed  pressure. 

The  vortex  is  spun  up  in  an  idealized  environment 
with  a  typical  tropical  sounding.  The  model  was  run  at 
a  horizontal  resolution  of  45  km,  and  on  an  /  plane. 
While  45 -km  horizontal  resolution  may  be  too  coarse  to 
resolve  the  inner-core  structure  of  some  TCs,  as  will  be 
shown,  a  reasonable  wind  structure  is  produced  in  com¬ 
parison  to  the  observations.  For  practical  application  of 
this  approach,  a  lookup  table  was  created  with  vortices 
that  were  spun  up  to  different  sea  level  pressures.  The 
environmental  pressure  was  set  to  a  typical  tropical  value 
of  1015  hPa  (Jordan  1958)  initially,  and  vortices  were 
created  with  targeted  surface  pressure  5-hPa  increments 
moving  to  900  hPa.  Experimentation  was  done  for  a  few 
cases,  varying  the  environmental  pressure  by  ±5  hPa  in 


Fig.  2.  Schematic  of  the  TCDI  system  flow.  Vortices  are  dy¬ 
namically  initialized  to  the  minimum  sea  level  pressure  (MSLP) 
using  the  idealized  model. 


the  ideal  model  to  account  for  variability  in  the  envi¬ 
ronmental  pressure  in  the  tropics.  It  was  found  that  it  did 
not  have  any  significant  effect  on  the  numerical  forecast 
in  COAMPS-TC  with  regard  to  intensity  and  track. 
However,  it  should  be  noted  that  large  variances  in  en¬ 
vironmental  pressure  would  affect  the  initial  vortex  in¬ 
tensity  significantly  because  the  integral  of  the  centrifugal 
and  Coriolis  terms  in  the  gradient  wind  equation  is  pro¬ 
portional  to  the  difference  between  the  environmental 
pressure  and  the  vortex  central  pressure.  We  also  wish  to 
note  that  the  ideal  vortex  is  blended  into  the  real  storm 
environmental  pressure  at  r  =  300  km  (see  next  para¬ 
graph),  and  the  environmental  pressure  in  the  ideal 
model  is  only  used  as  a  starting  point  to  spin  up  the 
vortex.  For  each  spunup  vortex,  a  48-h  time  integration  of 
the  model  was  used  to  allow  the  vortex  to  reach  an  ap¬ 
proximate  steady  state.  This  rather  short  time  period  for 
reaching  a  quasi-steady  state  could  be  used  since  each 
deeper  vortex  was  spun  up  beginning  with  the  vortex  that 
had  an  initial  sea  level  pressure  that  was  15  hPa  greater 
than  the  target  value. 

The  final  step  in  the  TCDI  procedure  is  to  insert  the 
new  vortex  back  into  the  filtered  model  initial  fields. 
This  was  accomplished  by  adding  the  new  vortex  per¬ 
turbations  (i.e.,  the  difference  of  the  total  field  with  the 
horizontal  mean  at  each  pressure  level)  to  the  filtered 
fields.  From  1013  to  500  mb,  a  blending  between  the 
TCDI  vortex  and  the  real  environment  is  accomplished 
by  taking  all  the  ideal  vortex  perturbations  for  r  <  200  km, 
followed  by  a  blending  zone  with  linear  weighting  from 
200  <  r  <  300  km,  and  taking  all  of  the  filtered  envi¬ 
ronmental  field  for  r  >  300  km.  Since  the  TC  upper- 
level  circulation  is  much  larger  than  at  lower  levels  (e.g., 
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Table  1.  North  Atlantic  (NA)  and  western  North  Pacific  (WNP)  tropical  cyclones  used  for  the  TCDI  testing. 


Basin 

Year 

TC 

Time  period 

NA 

2008 

Bertha  (02L) 

0000  UTC  6  Jul-1200  UTC  12  Jul 

NA 

2008 

Gustav  (07L) 

1200  UTC  25  Aug-0000  UTC  20  Sep 

NA 

2008 

Ike  (09L) 

1200  UTC  3  Sep-0000  UTC  9  Sep 

NA 

2009 

Bill  (03L) 

0000  UTC  16  Aug-0000  UTC  22  Aug 

NA 

2009 

Fred  (07L) 

0000  UTC  8  Sep-1200  UTC  11  Sep 

WNP 

2009 

Vamco  (11W) 

0000  UTC  18  Sep-1200  UTC  22  Aug 

WNP 

2009 

Choi- Wan  (15W) 

0000  UTC  13  Sep-1200  UTC  19  Sep 

WNP 

2009 

Parma  (19W) 

1200  UTC  29  Sep-0000  UTC  6  Oct 

WNP 

2009 

Melor  (20W) 

1200  UTC  29  Sep-0000  UTC  5  Oct 

WNP 

2009 

Lupit  (22W) 

1200  UTC  14  Oct-1200  UTC  24  Oct 

WNP 

2009 

Nida  (26W) 

1200  UTC  25  Nov-1200  UTC  29  Nov 

Elsberry  et  al.  1987),  this  blending  zone  slopes  linearly 
outward  with  decreasing  pressure  above  500  mb,  so  that 
at  p  =  10  hPa,  the  blending  zone  covers  the  region 
470  <  r  <  570  km.  When  running  a  mesoscale  model  with 
multiple  nests,  the  procedure  is  applied  to  each  nest,  with 
the  radial  blending  zone  increased  as  the  horizontal 


resolution  decreases  to  account  for  the  coarser  repre¬ 
sentation  of  the  TC  vortex  on  those  meshes. 

When  vortex  surgery  of  this  nature  is  done,  one  concern 
is  the  existence  of  gravity  wave  activity  caused  by  initial 
imbalances  in  the  model.  While  we  have  not  performed 
any  diagnostics  of  gravity  wave  activity,  it  is  possible  that 


a)  CONTROL 


b)TCDI 


Fig.  3.  Azimuthal  mean  structure  for  Bill  for  the  initial  conditions  of  the  0000  UTC  19  Aug  warm  start  using  the  (a)  control  and 
(b)  TCDI  methods,  (left)  The  azimuthal  mean  tangential  velocity,  (center)  the  azimuthal  mean  radial  velocity,  and  (right)  the  azimuthal 
mean  temperature  perturbation  (i.e.,  the  temperature  minus  the  area-average  temperature  on  each  pressure  level).  The  hypotenuse  of  the 
black  right  triangle  indicates  the  approximate  location  of  the  ocean  surface. 
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a)  CONTROL  b)TCDI 

t  =  0  h,  valid  0000  UTC  Aug  1 9  t  =  0  h,  valid  0000  UTC  Aug  1 9 


NHC  Best  Track:  955  hPa 

t  =  6  h,  valid  0600  UTC  Aug  1 9 


t  =  6  h,  valid  0600  UTC  Aug  1 9 


NHC  Best  Track:  954  hPa 


t  =  24  h,  valid  0000  UTC  Aug  20  t  =  24  h,  valid  0000  UTC  Aug  20 


Fig.  4.  Sea  level  pressure  evolution  (hPa)  for  the  0000  UTC  19  Aug  warm  start  of  Bill  for  the  (a)  control  and  (b)  TCDI  experiments  on  the 
5-km  COAMPS-TC  domain  with  a  horizontal  scale  of  approximately  905  km. 
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CONTROL^  =  0  h,  valid  0000  UTC  Aug  19 


CONTROL,  t=  24  h,  valid  0000  UTC  Aug  20 


TCDI, t  =  0  h,  valid  0000  UTC  Aug  1 9 


TCDI,t  =  24  h,  valid  0000  UTC  Aug  20 


Fig.  5.  Surface  wind  evolution  (kt)  for  the  0000  UTC  19  Aug  warm  start  of  Bill  for  (left)  the  control  and  (center)  TCDI  experiments  on 
the  5 -km  COAMPS-TC  domain  with  a  horizontal  scale  of  approximately  905  km.  (right)  The  H*Wind  surface  wind  analysis  is  shown  for 
approximately  the  same  time  (courtesy  of  NOAA/AOML/HRD). 


some  imbalances  do  exist  at  the  initial  time  causing  a  bal¬ 
ance  adjustment  and  associated  gravity  wave  activity.  In 
the  future,  it  is  planned  to  test  TCDI  with  a  digital  filter  to 
help  reduce  gravity  wave  activity  at  the  initial  time. 

Since  the  TCDI  procedure  only  modifies  the  vortex 
from  r  <  300  km  at  lower  to  middle  levels,  the  analysis  is 
retained  outside  of  this  radius.  Therefore,  the  initial  field 
with  the  original  synthetic  observations  at  4°  and  6°  radii 
is  still  included  with  the  TCDI  procedure,  yielding  an 
outer-wind  structure  that  closely  matches  the  observa¬ 
tions.  Thus,  TCDI  is  used  in  conjunction  with  3D  VAR  to 
help  improve  the  inner-core  intensity,  without  changing 
the  outer-wind  structure  from  the  analyses  generated  with 
the  TC  synthetic  observations. 

d.  Test  cases  and  best-track  data 

The  TCDI  method  was  tested  on  a  large  sample  of 
cases  from  the  2009  western  North  Pacific  and  2008/09 
North  Atlantic  seasons.  The  specific  TCs  and  date-time 
groups  used  to  test  the  TCDI  procedure  are  listed  in 


Table  1.  For  each  date-time  group,  two  72-h  forecasts 
were  executed:  the  control  forecast  and  the  forecast  using 
TCDI.  This  yielded  over  100  cases  for  most  lead  times. 

Intensity  and  track  comparisons  were  made  using  the 
NHC  best- track  data  in  the  North  Atlantic  basin  and  the 
Joint  Typhoon  Warning  Center  (JTWC)  best-track  data 
in  the  western  North  Pacific  basin. 

3.  Results 

In  this  section,  detailed  comparisons  of  the  TCDI  versus 
control  experiments  are  provided  for  two  cases  studies: 
North  Atlantic  Hurricane  Bill  (2009)  and  western  North 
Pacific  Typhoon  Choi- Wan  (2009).  Then,  the  intensity 
and  track  errors  are  compared  for  the  entire  sample. 

a.  Case  study:  Hurricane  Bill  (2009) 

Hurricane  Bill  was  the  most  intense  and  largest  tropical 
cyclone  of  the  relatively  inactive  2009  North  Atlantic 
hurricane  season.  It  formed  from  an  African  easterly  wave 
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a)  CONTROL 


b)TCDI 


BILL  Intensity:  2009081600 


BILL  Intensity:  2009081600 


BILL  Intensity:  2009081600 


BILL  Intensity:  2009081600 


Fig.  6.  Comparison  of  the  COAMPS-TC  intensity  forecasts  for  Hurricane  Bill  (2009).  (a)  The  control  run  intensity 
forecasts  and  (b)  the  TCDI  intensity  forecasts.  In  both  panels,  the  solid  black  line  is  the  NHC  best-track  intensity 
estimate  for  both  the  maximum  sustained  wind  and  MSLP,  and  individual  COAMPS-TC  72-h  forecasts  of  these 
quantities  are  shown  in  rainbow  colors.  Time  is  plotted  in  days  relative  to  0000  UTC  16  August. 


and  attained  tropical  depression  strength  at  0000  UTC 
15  August  2009.  Bill  intensified  steadily  and  become  a 
category  4  hurricane  with  maximum  sustained  winds  of 
115  kt  at  0600  UTC  19  August.  After  passing  approxi¬ 
mately  125  n  mi  west  of  Bermuda,  Bill  subsequently  re¬ 
curved  and  lost  its  tropical  characteristics  by  1200  UTC 
24  August  (Avila  2009). 

In  Fig.  3,  the  initial  conditions  for  the  0000  UTC  19 
August  warm  start  forecast  of  Bill  are  shown  for  both  the 
control  and  TCDI  methods.  At  this  time  Bill  had  NHC 
best-track  intensity  estimates  for  central  pressure  and  max¬ 
imum  sustained  winds  of  955  hPa  and  110  kt  (56  m  s_1), 
respectively.  As  shown  in  Fig.  3,  by  using  TCDI  a  higher 
azimuthal  mean  tangential  velocity  was  attained.  Ad¬ 
ditionally,  the  inflow  below  300  hPa  in  the  TCDI  initial 
conditions  is  more  realistic  than  the  weak  outflow  near  the 
core  region  in  the  control  experiment.  The  vortex  also  has 
a  more  pronounced  middle-  and  upper-level  warm-core 
structure  than  is  found  in  the  control  experiment,  and 
stronger  boundary  layer  inflow  and  upper-level  outflow. 


In  Fig.  4,  the  evolution  of  the  surface  pressure  field  is 
shown  for  the  control  and  TCDI  runs.  In  the  initial 
conditions,  the  control  experiment  central  pressure  is 
too  weak  (979  hPa),  while  the  TCDI  experiment  pro¬ 
duces  a  central  pressure  consistent  with  the  NHC  best- 
track  estimate  (955  hPa).  Also  note  that  the  horizontal 
structure  of  the  TCDI  run  is  realistic  and  matches  the 
control  experiment  in  the  outer  region.  During  the  time 
integration,  the  TCDI  run  stays  at  a  lower  central 
pressure  (closer  to  the  NHC  best-track  intensity),  while 
the  control  run  remains  too  weak.  In  Fig.  5,  a  compari¬ 
son  of  the  10-m  wind  field  evolution  is  shown  for  the 
control  and  TCDI  runs  at  the  initial  time  and  at  t  =  24  h. 
Also  shown  in  the  right  panels  of  Fig.  5  is  the  H*Wind 
surface  wind  analysis  (Powell  and  Houston  1996;  Powell 
and  Houston  1998;  Powell  et  al.  1998)  at  approximately 
the  same  time  [courtesy  of  the  National  Oceanic  and 
Atmospheric  Administration/ Atlantic  Oceanographic  and 
Meteorological  Laboratory/Hurricane  Research  Division 
(NOAA/AOML/HRD)].  At  the  initial  time,  qualitatively 
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^CONTROL 
t  =  24  h 


Fig.  7.  Comparison  of  the  observed  structure  of  Supertyphoon  Choi-wan  with  the  control  and  TCDI  numerical  experiments  with 
COAMPS-TC.  The  COAMPS-TC  24-h  forecast  is  valid  at  0000  UTC 14  Sep  and  the  48-h  forecast  is  valid  at  0000  UTC  15  Sep.  (left)  The 
control-experiment-simulated  radar  reflectivity  at  t  =  0  h  and  t  =  48  h,  (middle)  the  TCDI  experiment  simulated  radar  reflectivity  at  t  =  0  h 
and  t  =  48  h,  and  (right)  the  85-/89-GHz  microwave  data,  courtesy  of  the  Naval  Research  Laboratory. 


both  the  control  and  TCDI  experiments  have  a  structure 
consistent  with  the  H*Wind  analysis  with  regard  to  the 
radii  of  50-  and  34-kt  winds.  The  radius  of  maximum  wind 
is  slightly  larger  in  both  the  control  and  TCDI  runs 
than  in  the  H*Wind  analysis.  At  t  =  24  h,  both  the 
control  and  TCDI  produce  reasonable  structures  in 
comparison  to  the  H*Wind  analysis  (including  the  wave- 
number  one  asymmetry  due  to  the  storm  motion  and  the 
radii  of  50-  and  34-kt  winds).  At  this  time,  the  TCDI 
run  produces  a  more  intense  vortex  consistent  with  the 
H*Wind  analysis. 

The  intensity  forecast  comparison  between  TCDI  and 
the  control  method  for  all  of  the  COAMPS-TC  forecasts 
for  Hurricane  Bill  is  shown  in  Fig.  6.  In  the  left  panels  in 
Fig.  6,  the  performance  of  the  control  experiment  is 
shown,  and  in  the  right  panels  the  performance  of  the 
TCDI  experiment  is  shown.  In  both  sets  of  panels  the  solid 
black  lines  are  the  NHC  best-track  estimates  of  the  min¬ 
imum  central  pressure  and  maximum  sustained  surface 
wind.  Individual  72-h  COAMPS-TC  forecasts  of  these 
quantities  are  depicted  as  colored  curves,  with  the  first 
cold-start  forecast  occurring  at  0000  UTC  16  August,  and 
subsequent  warm  starts  every  12  h  thereafter.  The  ordi¬ 
nate  shows  days  relative  to  0000  UTC  16  August. 


In  comparing  the  pressure  plots,  it  is  readily  apparent 
that  the  3D  VAR  system  with  synthetic  observations  in 
the  control  experiment  is  generally  not  able  to  initialize 
Bill  at  the  observed  intensity,  being  too  weak  by  both 
parameters.  In  contrast,  by  using  TCDI,  the  model  is 
initialized  very  close  to  the  observed  surface  pressure. 
The  improvement  in  the  initial  condition  intensity  has 
a  lasting  effect  on  the  subsequent  evolution  of  the  vor¬ 
tex.  By  initializing  at  the  correct  intensity,  the  intensity 
forecast  using  TCDI  is  significantly  better  than  that  of 
the  control  experiment.  A  similar  result  is  seen  in  the 
maximum  wind  plots:  the  control  experiment  is  gener¬ 
ally  initialized  too  weak,  while  the  TCDI  initial  intensity 
is  closer  to  the  observations.  As  shown  in  Fig.  6,  the  in¬ 
tensity  trend  is  generally  similar  between  the  control  and 
TCDI  experiments,  and  the  intensity  error  reduction  in 
the  TCDI  experiments  is  largely  due  to  the  TC  being  ini¬ 
tialized  very  close  to  the  NHC  best-track  intensity.  Both 
the  control  and  TCDI  runs  suffer  from  a  spindown  in  the 
maximum  sustained  wind  after  initialization.  Since  TCs 
typically  have  a  horizontal  scale  smaller  than  the  vortex 
Rossby  length,  if  there  is  an  initial  imbalance,  the  mass 
field  will  typically  adjust  to  the  wind  field.  Therefore,  it  is 
not  clear  that  the  spindown  is  directly  related  to  balance 
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a)  CONTROL  b)  TCDI 

CHOI-WAN  Intensity:  2009091112  CHOI-WAN  Intensity:  2009091112 


Fig.  8.  Comparison  of  the  COAMPS-TC  intensity  forecasts  for  Supertyphoon  Choi- wan  (2009).  (a)  The  control 
run  intensity  forecasts  and  (b)  the  TCDI  intensity  forecasts.  In  both  panels,  the  solid  black  line  is  the  NHC  best-track 
intensity  estimate  for  both  the  maximum  sustained  wind  and  MSLP,  and  individual  COAMPS-TC  72-h  forecasts  of 
these  quantities  are  shown  in  rainbow  colors.  Time  is  plotted  in  days  relative  to  1200  UTC  11  September. 


issues.  One  possibility  is  that  the  analysis  increments  for 
both  the  control  and  TCDI  are  too  large  for  strong  storms, 
resulting  in  a  tendency  to  move  toward  the  model  attractor 
(i.e.,  the  first-guess  field)  after  initialization.  It  is  also  pos¬ 
sible  that  boundary  layer  processes  are  causing  the  surface 
wind  to  spin  down  initially.  This  topic  will  be  investigated 
more  rigorously  in  the  future. 

b.  Case  study:  Supertyphoon  Choi-wan  (2009) 

A  tropical  depression  formed  from  an  area  of  con- 
vectively  unstable  environment  in  the  western  North 
Pacific  Ocean  approximately  1100  km  west  of  Guam  on 
1800  UTC  11  September  2009.  Shortly  thereafter,  the  de¬ 
pression  strengthened  to  a  tropical  storm  and  was  named 
Choi-wan.  Choi-wan  subsequently  underwent  a  number 
of  rapid  intensification  events  and  was  classified  as  a 
supertyphoon  by  JTWC  at  1800  UTC  14  September  with 
maximum  sustained  surface  winds  of  130  kt.  Choi-wan 
intensified  further  to  140  kt  and  918  hPa  at  1800  UTC  15 
September  and  remained  a  supertyphoon  until  0000  UTC 


17  September.  Afterward,  Choi-wan  underwent  an  eye- 
wall  replacement  cycle,  resulting  in  minor  weakening. 
More  significant  weakening  followed  as  it  entered  an  en¬ 
vironment  of  strong  shear.  Following  that,  it  recurved 
and  transitioned  into  an  extratropical  cyclone. 

In  Fig.  7,  horizontal  plots  of  the  control  and  TCDI 
experiment  radar  reflectivities  are  shown  at  t  =  24  and 
48  h  from  the  cold  start  initialized  at  0000  UTC  13 
September,  along  with  the  Tropical  Rainfall  Measuring 
Mission  (TRMM)  Microwave  Imager  (TMI)  85-GHz 
channel  and  Advanced  Microwave  Scanning  Radiometer 
(AMSR-E)  89-GHz  channel  brightness  temperatures  in 
the  vicinity  of  Choi- wan  at  times  closely  corresponding  to 
the  forecast  times.  The  horizontal  scale  of  all  plots  was 
made  to  be  exactly  the  same  to  accurately  compare  their 
structures.  In  the  TMI- AMSR-E  plots,  brighter  colors 
denote  a  reduction  in  brightness  temperature  due  to 
more  scattering  by  ice  hydrometeors  in  deep  convection, 
and  can  be  qualitatively  compared  to  the  COAMPS-TC- 
simulated  radar  reflectivity. 
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Maximum  Sustained  Wind  Statistics  (Warm  Starts) 


Minimum  Central  Pressure  Statistics  (Warm  Starts) 


N  =  (112)  (112)  (111)  (109)  (106)  (104)  (97) 


Fig.  9.  Intensity  mean  absolute  errors  and  biases  by  (top)  max¬ 
imum  sustained  wind  and  (bottom)  minimum  central  pressure  for 
the  homogenous  sample  of  warm-start  cases.  The  number  of  cases 
at  each  lead  time  (in  intervals  of  12  h)  is  shown  in  parentheses.  The 
thin  black  line  indicates  zero. 

Neither  the  control  nor  the  TCDI  run  is  able  to  cap¬ 
ture  the  intricate  details  of  the  observed  inner  core  at 
48  h.  Such  intricate  details  would  be  very  difficult  to 
predict  in  a  mesoscale  model  even  at  very  high  horizontal 
resolution.  These  experiments  (at  5-km  horizontal  reso¬ 
lution)  are  certainly  too  coarse  to  adequately  capture 
these  details.  However,  at  both  lead  times,  the  structure 
of  the  TCDI  run  is  more  consistent  with  the  observations 
than  the  control  run.  The  TCDI  run  produces  a  stronger 
vortex,  which  more  closely  matches  the  observed  track 
and  intensity  (not  shown),  and  this  stronger  vortex 
has  an  inner-core  structure  that  more  closely  matches  the 
observations.  Note  in  particular  that  the  TCDI  run  at 
t  =  24  h  has  a  more  distinct  eye  and  eyewall,  and  an  az¬ 
imuthal  wavenumber  one  asymmetry,  which  was  ob¬ 
served  from  TMI  (although  approximately  90°  out  of 
phase).  At  t  =  48  h,  the  TCDI  run  produces  a  tighter 
eyewall  than  the  control  run,  which  is  more  consistent 
with  the  TMI  observations. 


Fig.  10.  Track  errors  for  warm-start  cases  using  the  control  and 
TCDI  methods.  The  number  of  cases  at  each  lead  time  (in  intervals 
of  12  h)  is  shown  in  parentheses. 


The  intensity  forecast  comparison  between  TCDI  and 
the  control  method  for  all  of  the  COAMPS-TC  forecasts 
for  Supertyphoon  Choi-wan  is  shown  in  Fig.  8.  In  the  left 
panels  in  Fig.  8,  the  performance  of  the  control  experi¬ 
ment  is  shown,  and  in  the  right  panels  the  performance 
of  the  TCDI  experiment  is  shown.  In  both  sets  of  panels 
the  solid  black  lines  are  the  JTWC  best-track  estimates 
of  the  minimum  central  pressure  and  maximum  sustained 
surface  wind.  Individual  72-h  COAMPS-TC  forecasts  of 
these  quantities  are  depicted  as  colored  curves,  with  the 
first  cold-start  forecast  occurring  at  1200  UTC  11  Sep¬ 
tember,  and  subsequent  warm  starts  every  12  h  there¬ 
after.  The  ordinate  shows  days  relative  to  1200  UTC  11 
September.  In  the  pressure  plots  (Fig.  8,  top),  the  control 
cases  of  Choi-wan  are  generally  initialized  too  weak  in 
comparison  to  the  JTWC  best-track  estimate.  However, 
with  the  TCDI  method,  the  TC  is  initialized  much  deeper, 
closer  to  the  best-track  estimate.  For  very  deep  pressures 
(<930  hPa),  the  TCDI  run  fills  initially  and  then  deepens 
afterward.  This  initial  filling  may  be  due  to  the  model 
attractor  issue  discussed  previously  for  Hurricane  Bill. 
By  maximum  sustained  wind  (Fig.  8,  bottom),  the  TCDI 
runs  generally  perform  better  than  the  control  runs. 
Typically,  the  TC  is  initialized  closer  to  the  JTWC  best- 
track  estimate,  and  there  is  less  of  a  spindown  in  TCDI 
as  compared  to  the  control  runs.  In  summary,  similar  to 
Hurricane  Bill,  the  TCDI  run  exhibited  a  more  accurate 
initial  intensity  and  intensity  forecast  than  the  control 
run  as  compared  to  the  JTWC  best-track  data  for  most 
of  the  forecasts. 

c.  Homogenous  sample  statistics 

In  Fig.  9,  the  intensity  statistics  for  the  homogenous 
sample  are  shown  for  all  warm-start  cases,  as  listed  in 
Table  1.  The  number  of  cases  is  listed  in  parentheses  at 
each  lead  time  from  t  =  0, 12,  24, . . .  ,  to  72  h.  The  black 
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Maximum  Sustained  Wind  Statistics  (Cold  Starts) 


Minimum  Central  Pressure  Statistics  (Cold  Starts) 


Fig.  11.  Intensity  mean  absolute  errors  and  biases  by  (top) 
maximum  sustained  wind  and  (bottom)  minimum  central  pressure 
for  the  cold-start  cases.  The  number  of  cases  at  each  lead  time  (in 
intervals  of  12  h)  is  shown  in  parentheses.  The  thin  black  line  in¬ 
dicates  zero. 

curves  depict  the  control  experiments  and  the  gray 
curves  depict  the  TCDI  experiments.  The  solid  lines  are 
the  mean  absolute  errors  and  the  dashed  lines  are  the 
biases.  For  all  lead  times,  the  TCDI  experiments  have 
mean  absolute  errors  in  maximum  sustained  surface  wind 
that  are  approximately  5  kt  lower  than  the  control  ex¬ 
periments.  Although  both  datasets  have  a  negative  bias 
(i.e.,  the  intensity  is  generally  too  weak),  this  bias  is  re¬ 
duced  in  the  TCDI  runs.  This  is  due  to  the  fact  that  in 
general  the  TCDI  runs  are  initialized  closer  to  the  best- 
track  intensity  and  stay  more  intense  during  the  forecast. 
A  similar  result  is  seen  by  the  minimum  central  pressure 
statistics  (Fig.  9,  bottom).  The  TCDI  run  mean  absolute 
error  and  bias  are  greatly  reduced  at  t  =  0  h,  which  is 
expected  since  the  TCDI  method  uses  nudging  to  the 
observed  surface  pressure.  At  subsequent  lead  times, 
the  TCDI  experiments  show  reduced  errors  and  biases 
with  regard  to  surface  pressure  as  well.  The  track  error 
comparison  of  both  methods  for  warm-start  cases  is 


Fig.  12.  Track  errors  for  cold-start  cases  using  the  control  and 
TCDI  methods.  The  number  of  cases  at  each  lead  time  (in  intervals 
of  12  h)  is  shown  in  parentheses. 


shown  in  Fig.  10.  There  were  found  to  be  no  significant 
differences  in  track  error  between  the  two  groups. 

In  Fig.  11,  the  intensity  statistics  for  all  cold-start  cases 
as  listed  in  Table  1  are  depicted.  Since  there  are  11  TCs, 
there  are  11  cold  starts  in  which  the  first  TC  forecast  uses 
the  NOGAPS  global  analysis  for  the  first  guess.  By 
maximum  sustained  wind,  the  TCDI  method  reduces  the 
bias;  however,  no  improvement  was  seen  in  the  error.  In 
contrast,  significant  improvements  were  seen  in  the  error 
and  bias  with  regard  to  the  minimum  central  pressure  by 
using  TCDI.  There  were  also  no  significant  differences 
in  track  errors  by  using  either  method  (Fig.  12).  The 
TCDI  method  performed  slightly  worse  initially,  and 
slightly  better  between  48  and  72  h;  however,  these 
differences  are  not  statistically  significant.  The  lack  of 
improvement  in  the  intensity  error  by  maximum  sus¬ 
tained  wind  for  cold-start  cases  is  somewhat  surprising. 
However,  since  the  current  TCDI  method  uses  sea  level 
pressure  nudging,  the  TCDI  wind-pressure  relationship 
may  differ  slightly  from  the  observed  TC  wind-pressure 
relationship.  Since  the  wind  bias  is  negative,  the  TC  size 
using  TCDI  is  likely  smaller  than  the  observations  for 
these  11  cases.  It  should  be  noted,  however,  that  11  cases 
is  not  a  large  enough  sample  for  statistical  significance. 

4.  Summary 

A  tropical  cyclone  dynamical  initialization  (TCDI) 
method  applicable  to  TC  numerical  prediction  systems 
is  described  and  its  performance  is  tested  on  over  100 
cases  in  2008  and  2009  in  the  North  Atlantic  and  western 
North  Pacific  Ocean  basins  using  the  high-resolution 
full-physics  COAMPS-TC  model.  The  method  involves 
the  removal  of  the  analyzed  TC  vortex  and  the  insertion 
of  a  new  vortex  that  is  dynamically  spunup  to  the  ob¬ 
served  surface  pressure  using  an  independent  full-physics 
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primitive  equation  model.  The  results  indicate  that  use 
of  such  a  method  has  potential  for  both  a  reduction  in 
intensity  errors  and  biases,  and  improved  initial  and 
forecasted  TC  structures,  with  no  degradation  in  track 
performance.  For  the  sample  of  cases  here,  maximum 
1-min  sustained  wind  intensity  errors  were  reduced  by  ap¬ 
proximately  5  kt  over  all  lead  times  from  t  =  0, 12, 24, ... , 
to  72  h  by  using  TCDI.  The  primary  reason  for  the  re¬ 
duction  in  intensity  errors  is  due  to  the  TC  vortex  being 
initialized  close  to  the  best-track  intensity.  In  addition, 
the  dynamic  initialization  method  also  produces  a  self- 
consistent,  balanced  vortex,  with  realistic  boundary  and 
outflow  layers,  potentially  helping  to  improve  the  forecast. 

The  dynamic  initialization  procedure  has  long  been 
developed  and  used  in  the  GFDL  model.  One  major 
difference  between  the  GFDL  system  and  our  system  is 
that  GFDL  uses  an  axisymmetric  version  of  its  model  for 
the  spinup  of  the  vortex  while  we  use  a  different,  inde¬ 
pendent,  three-dimensional  primitive  equation  model  for 
the  spinup.  Therefore,  we  demonstrate  that  it  is  not  nec¬ 
essary  to  use  the  same  model  as  the  forecast  model  for 
spinup,  as  the  goal  is  to  have  a  TC  vortex  that  is  more 
complete  three-dimensionally  and  more  balanced  in  the 
vertical  with  realistic  outflow  at  upper  levels  and  boundary 
layer  inflow.  Our  point  is  that,  for  any  TC  prediction  sys¬ 
tem,  one  can  take  the  TCM3  or  GFDL  idealized  version 
and  use  it  for  TCDI  to  improve  the  initial  state  of  a  storm. 

There  are  a  number  of  future  improvements  that  can 
be  incorporated  into  the  TCDI  system.  First,  the  im¬ 
plementation  of  wind  structure  dynamic  initialization 
will  allow  for  the  relaxation  to  a  prescribed  radial  wind 
profile,  such  as  a  modified  Rankine  vortex.  This  will 
allow  initialization  of  TCs  with  the  correct  inner-core 
structure  as  well  as  intensity.  Second,  the  use  of  TRMM 
satellite-derived  diabatic  heating  profiles  as  forcing  could 
further  enhance  the  ability  to  more  closely  capture  the 
observed  TC  asymmetric  structure.  Finally,  while  we 
have  shown  this  system  to  be  portable  in  that  the  pre¬ 
diction  model  and  idealized  TCDI  model  can  be  dif¬ 
ferent,  in  the  future  it  is  planned  to  use  the  idealized 
version  of  COAMPS-TC  for  TCDI.  This  may  yield  a  slight 
improvement  in  performance  since  both  models  are  based 
on  the  same  governing  equations,  physics,  and  numerics. 

In  recent  workshops  on  NOAA‘s  Hurricane  Forecast 
Improvement  Project  (HFIP),  the  proper  initialization 
of  TCs  in  numerical  prediction  systems  remains  one  of 
the  most  challenging  tasks  for  TC  intensity  prediction. 
TCDI,  as  envisioned  here,  can  be  a  temporary  improve¬ 
ment  that  can  be  included  in  a  3D  VAR  initialization 
system  until  4DVAR,  ensemble  Kalman  filter,  and  hy¬ 
brid  mesoscale  data  assimilation  techniques  in  opera¬ 
tional  numerical  TC  prediction  systems  become  more 
feasible. 
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